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The interaction of ultra-intense lasers with solid foils can be used to accelerate ions to high 
energies well exceeding 60 MeV [1]. The non-linear relativistic motion of electrons in the intense 
laser radiation leads to their acceleration and later to the acceleration of ions. Ions can be accelerated 
from the front surface, the foil interior region, and the foil rear surface (TNSA, most widely used), 
or the foil may be accelerated as a whole if sufficiently thin (RPA). 

Here, we focus on the most widely used mechanism for laser ion-acceleration of TNSA. Starting 
from perfectly flat foils we show by simulations how electron filamentation at or inside the solid 
leads to spatial modulations in the ions. The exact dynamics depend very sensitively on the chosen 
initial parameters which has a tremendous effect on electron dynamics. In the case of step-like 
density gradients we find evidence that suggests a two-surface-plasmon decay of plasma oscillations 
triggering a Raileigh-Taylor-like instability. 


In laser ion acceleration, an ultra intense laser inter¬ 
acts with a solid and quickly ionizes its surface. The laser 
light creates an immense pressure on the surface, pro¬ 
duces large and fast electron currents into the solid and 
excites plasma waves. In such a scenario many instabili¬ 
ties develop which can disrupt the surface and break up 
the electron currents. 

It has been shown previously that instabilities in ultra- 
thin solid foils can imprint onto the spatial structure of 
ions in the regime of relativistically induced transparency 
target normal sheath acceleration (RIT TNSA |2|) or ra¬ 
diation pressure acceleration (RPA [3]) [4] ~ mostly at¬ 
tributed to the Rayleigh-Taylor instability (RTI) [[5]. In¬ 
stabilities behind the foils can also develop during the ps 
timescales of the plasma expansion of TNSA [6]. Here we 
show that instabilities developing in solid foils can also 
be of significance for the laser absorption and ion accel¬ 
eration for thicker foils where no RIT occurs. We present 
results from particle-in-cell (PIC) simulations using the 
code iPlCLS2D and refer the reader to our Ref. [8] 
for experimental evidence. We use dimensionless units 
with c = e = m e = cjo = 1 and measure distances in 
x-direction (the laser direction) relative to the initial foil 
front surface position and times relative to the time when 
the laser intensity maximum reaches x = 0. In most of 
our analysis we will focus on two simulations, compar¬ 
ing the cases of a flat foil with exponential preformed 
plasma on the surface (simulation A) and without (simu¬ 
lation B), the full simulation parameters are given in the 
caption of Fig. [lj 

One main finding is that instabilities which are prone 
to evolve in all of our simulations imprint on the ion 
density distribution behind the foil slab even for the mi¬ 
cron thick foils used. The instabilities we discuss in this 
paper are developing in the electrons first since the heav¬ 


ier ions remain almost at rest during the relevant time 
scales of electron motion. Of course this is not strictly 
correct since otherwise no RTI instability could develop. 
Fig. QJa) shows exemplary part of the electron density 
distribution of simulation At = 187r = 9 To after the 
laser maximum has reached the foil front surface, and 
Fig. QJb) shows the corresponding ion density distribu¬ 
tion. In this simulation the electrons pushed into the foil 
at 2 cjo by the laser Lorentz force are prone to a transverse 
Weibel-like instability Q. The 2ujq electron slabs become 
filamented transversely and create energetic channels al¬ 
most normal to the foil surface which are surrounded by 
strong quasi-static magnetic fields. This happens during 
the first few plasma wavelengths X p = 2ttuj~ 1 = 2im e ^ 2 
as can be seen in Fig. [2^. The ion acceleration at the foil 
rear surface then follows this structure(cp. Fig, []Ji,b). 
It is this that is different to the mechanism described 
in [6]: the rear surface ion’s acceleration itself is struc¬ 
tured transversely from the beginning since the driving 
energetic electrons reach the rear surface filamented al¬ 
ready in contrast to the relatively slow Weibel instability 
witnessed in [6] during the ion acceleration. Of course 
both can happen concurrently. 

We now study the case of step-like density gradient at 
the front surface (simulation B). In comparison to simula¬ 
tion A, where the bulk electrons remained distributed al¬ 
most homogeneously transversely at the surface and only 
energetic electrons became filamented inside the foil, here 
the foil surface already shows transverse structure both in 
electron density (Fig. [2Jd,c) and energy density (Fig.[2k), 
with the latter not increasing during passage through the 
foil contrarily to simulation A. The reason for this lies in 
the density ripples which give rise to filamented injection 
of electrons, similar to the mechanism described in 0- 
We therefore conclude that in simulation B another fila- 
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FIG. 1. (a) Part of the electron energy density distribution at t — 18tt = 9 To after the laser , normalized to 2/3 of its maximum 
and (b) ion density distribution of simulation A, and (c) electron energy distribution of simulation B. Simulation parameters: 
only the plane defined by the laser propagation direction (x-axis) and the electric held polarization (y-axis) is simulated, 
assuming invariance of the system in z-direction; simulation box 10Ao x 10Ao, 192 cells per Ao and laser period To; laser is 
modeled by a transversely plane wave with short gaussian rise with 1 sigma-width of 2Ao/c followed by a hat top, wavelength Ao, 
peak intensity Iq [W/cm 2 ] = 1.38 X 10 2O /(Ao [(im]) 2 corresponding to a normalized held strength ao = [2Io/(n c m e c 3 )] 1//2 = 10 
where n c = m. e eouJoe 2 is the critical electron density for the laser with angular frequency o;o; laser is aligned normal to the 
target foil surface; hat target foils modeled by a fully preionized plasma slab of 2Ao thickness: ions with charge-to-mass ratio 
Q/M — 1/2 and neutralizing electrons with density n e ,o = 100 n c ; front of the foils is modeled by an exponentially increasing 
preplasma with scale length L, (a) and (b) L — 27r/10 = 0.1 Ao (simulation A) and (c) L — 0 (simulation B). 



FIG. 2. (a) Transverse electron energy density contrast when co-moving with one of the 2ujq electron slabs after it is injected 
into the solid for simulation A (black line) and simulation B (gray/red line), (b, c) Electron density snapshot at t — 18 tv = 9 To 
for (b) simulation A and (c) simulation B. 


mentation mechanism must be present. Density ripples 
usually are created by RTI and the density distribution in 
Fig. [2t shows a pattern consistent with the linear stage of 
RTI. Yet, typically the average distance A s between den¬ 
sity maxima (wave number ks = A^ 1 ) should be given 
by the smallest transverse scale that could seed the in¬ 
stability, since the growth rate T = (gks) -1 ^ 2 (g is the 
acceleration of the surface) is largest for the largest ks- 
Hence we are looking for a mechanism that can seed a 
transverse density fluctuation which is then prone to a 
RTI growth |111, [l2|. Such a mechanism was described 
e.g. in 0 [ 15 ] where a generation of 2D electron surface 
oscillations (2DESOs) by a parametric decay of 2ujq os¬ 
cillations of the foil surface in the laser field was found. 
There from initially purely longitudinal driving oscilla¬ 
tions of the surface (for the 2co>o oscillations the driver 
frequency is ud = 2cjo), oscillations in a standing wave 
along the surface are generated via excitation of upward 
and downward moving and interfering surface waves with 
us = ud /2 and ks given by [Eqn. (6) in 0 ] . For our 
simulation conditions this would correspond to ks — 1, 
much smaller than the value extracted from the PIC 
simulation k$ IC = 7.7. This decay of 2ujq oscillations 


with the characteristic Icjo longitudinal oscillation is ob¬ 
served in our simulation only at later times growing only 
slowly, which we attribute to larger ks- However, con¬ 
sidering not the 2 cjo oscillation as a driving source but 
the plasma oscillations excited by the disturbances we 
find good numerical agreement. Fig. [3] shows a compar¬ 
ison of kg IC for a number of simulations we performed 
with different parameters (varying n e? o and ao) with the 
analytical value ks- The agreement is generally rather 
good. Also, we find oscillations at uo p /2 in the spectrum 
of surface oscillations, Fig. [3^b). Here we defined a den¬ 
sity level and recorded the position in x-direction (laser 
direction) where the rising density reaches that value. 
We then computed the temporal Fourier transform of 
this position averaged over the positions along the y- 
direction at the density ripple maxima, neglecting the 
regions in between. The result would be a strong signal 
around the laser harmonics (overlying all other signals) 
which however are expected to be transversely flat here. 
Therefore we took the oscillation at the ” ‘inner surface” ’. 
By this we mean the oscillation of the falling density 
slope behind the region of laser pressure steepened elec¬ 
tron density, which can screen the lower laser harmonics. 
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FIG. 3. (a) kg IC versus calculated ks from [Eqn. (6) in [l5[] not including thermal effects, (b) and (c) Temporal Fourier 
transform of surface oscillations for (b) simulation B and (c) a simulation with same parameters as B but n e ,o = 300 and 
clq = 30 (details see text). The time over which the oscillations were evaluated is indicted in the legend. 



FIG. 4. Energy density distribution for simulations with pa¬ 
rameters as simulation B but with initial surface roughness. 
The largest spatial frequency is indicated by the black bar, the 
smallest is twice this size. The amplitude of the surface rough¬ 
ness is O.lAo. Color scale is the same as in Fig. [I] Largest 
surface ripple growth occurs when the spatial frequency ks 
is seeded, with a characteristic mushroom shape (non-linear 
RTI) growing at the front (central panel). Here, kg 10 = 7.6, 
k s = 7.3. 


As expected we find a signal at the plasma frequency 
uj p = (n e ,m ax /Tmax) 1 / 2 computed using the maximum 
density and temperature at the peak of the laser pres¬ 
sure steepened density profile averaged over the distance 
of a plasma period. Additionally we also found another 
pronounced peak at uj p /2. Since there is still a significant 
contribution from laser harmonics, we performed another 
simulation with n e? o = 300 and ao = 30 (Fig. [3{c)). Due 
to the higher density the skin depth S = cuj~ x is much 
shorter, screening laser harmonics, and additionally the 
plasma frequency and its half value are shifted to higher 
frequencies and hence are at higher and hence weaker 
harmonics. Here we again obtain a very clear signal at 
uj p and ujp/2. 

It is interesting to mention that the development of sur¬ 
face ripples seen in simulation B can be artificially seeded. 
Introducing an initial surface roughness containing a mix¬ 
ture of several spatial frequencies around ks = 2 ksw we 
first find that the smallest spatial frequencies have the 
fastest growth, as expected for RTI, and that this growth 
is most dominant when the spatial frequencies are initial¬ 
ized around ks (Fig.0J). For yet higher spatial frequencies 
the effect of seeding abruptly ceases. 

Summarizing, we found by simulations that the laser 
ion acceleration process in the TNSA regime can be in¬ 
fluenced by electron instabilities at the surface or inside 


the foil. The exact mechanism and scale of spatial mod¬ 
ulations in the beam of hot electrons and ions sensitively 
depends on the initial parameters. For the pursuit of 
higher ion energies at higher laser intensities, e.g. avail¬ 
able at (future) Petawatt laser systems, a deeper under¬ 
standing and further characterization of laser and target 
parameter dependent plasma instabilities will be essen¬ 
tial. 
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